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Abstract 

We introduce a class of singular partial differential equations, the second-order hyperbolic 
Fuchsian systems, and we investigate the associated initial value problem when data are 
imposed on the singularity. First of all, we analyze a class of equations in which hyperbolicity 
is not assumed and we construct asymptotic solutions of arbitrary order. Second, for the 
proposed class of second-order hyperbolic Fuchsian systems, we establish the existence of 
solutions with prescribed asymptotic behavior on the singularity. Our proof is based on 
a new scheme which is also suitable to design numerical approximations. Furthermore, as 
shown in a follow-up paper, the second-order Fuchsian framework is appropriate to handle 
Einstein's field equations for Gowdy symmetric spacetimes and allows us to recover (and 
slightly generalize) earlier results by Rendall and collaborators, while providing a direct 
approach leading to accurate numerical solutions. The proposed framework is also robust 
enough to encompass matter models arising in general relativity. 
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1 Introduction 

This is the first of a series of papers [7j [8] devoted to solving the initial value problem for certain 
classes of spacetimes of general relativity. Specifically, we are interested in spacetimes enjoying 
certain symmetries, especially the Gowdy symmetry, and in a formulation when data are imposed 
on a singular hypersurface where curvature generically blows-up. 

For instance, one may consider (3 + l)-dimensional, vacuum spacetimes (M,g) with spatial 
topology T 3 , satisfying the vacuum Einstein equations under the Gowdy symmetry assumption, 
i.e. under the existence of an Abelian T 2 isometry group with spacelike orbits whose so-called 
"twist constants" [5] vanish identically. These so-called Gowdy spacetimes on T 3 were first 
studied in |12| . A combination of theoretical and numerical works has led to a detailed picture 
of the behavior of solutions to the Einstein equations as one approaches the singular boundary 
of such spacetimes; see [5| [13 ] US ) l2" H W2 \ l2"5] . 

For the above analysis, one important tool was provided by Rendall and his collaborators 
[UITJIH!] who developed the so-called Fuchsian method to handle the singular evolution equations 
associated with the Einstein equations for such spacetimes. This method allowed the authors 
to derive precise information about the behavior of solutions near the singularity, which was a 
key step in the general proof of Penrose's strong cosmic conjecture eventually established by 
Ringstrom [23] . 

Our aim here and in the follow-up papers [7JIB] is two-fold. On one hand, we re- visit Rendall's 
theory (which covers smooth solutions to first-order equations) and we develop here a well- 
posedness theory in Sobolev spaces for the class of second-order hyperbolic Fuchsian systems, 
defined below in Section 3. As we will show, this class includes many systems of equations arising 
in general relativity. 

On the other hand, following a strategy initially proposed in Amorim, Bcrnardi, and LeFloch 
[T] , we are interested in the numerical approximation of these Fuchsian equations when data are 
imposed on the singularity of the spacetime and one evolves the solution from the singularity. In 
contrast, standard numerical approaches consider the evolution toward the singularity. Our main 
improvement here upon [T] is that, in short, no restriction need be imposed on the coefficients of 
the Fuchsian system, provided asymptotic expansions of arbitrary large order are sought. This 
issue will be developed in [TJ. 

The present paper is theoretical in nature, and our main contributions can be summarized 
as follows: 

• Second-order formulation. First-order Fuchsian systems have been used successfully 
to handle the equations describing (vacuum) Gowdy spacetimes [H [TJl [15]. However, 
we argue here that second-order hyperbolic Fuchsian systems, as we define them in this 
paper, arise more naturally in the applications. For instance, Gowdy spacetimes described 
by second-order partial differential equations, and it is natural to keep the second-order 
structure. In particular, expansions required in the theoretical analysis (as well as in 
actual computations required for the numerical discretization) are also more natural with 
the second-order formulation. 

• Hyperbolicity property. In addition, in the applications to general relativity, ensuring 
and checking the hyperbolicity of the equations under consideration (after a suitable re- 
duction of the original equations) is expected to be more convenient with the second-order 
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formulation. For instance, it is easier to recognize from the original second-order system if 
the equations form a system of coupled wave equations, while this property is much less 
evident in the first-order formulation. As the discussion in [TO] demonstrates, it can be 
cumbersome to formulate a system of coupled non-linear wave equations as a first-order 
hyperbolic Fuchsian system satisfying all the properties required for local well-posedness. 

• Singular part. We also introduce here a construction algorithm which includes the singu- 
lar part of the solution. This is different from the classical first-order Fuchsian approach, 
where one first makes an ansatz by removing from the solution its (expected) singular part, 
and then applies the Fuchsian theory to the (regular) remainder. Our approach leads to the 
notion of the singular initial value problem for Fuchsian hyperbolic equations which can be 
understood as a generalization of the initial value problem for (standard) hyperbolic equa- 
tions. This singular initial value problem covers cases whose singularity is oscillatory in 
nature, in a manner consistent with the so-called BKL conjecture (introduced by Bclinsky, 
Khalatnikov, and Lifshitz). 

We shall, first, define the class of systems of interest, that is, the second-order (hyperbolic) 
Fuchsian systems and, then, investigate the associated initial value problem when data are im- 
posed on the singularity. Precisely, in Section 2, we analyze a class of equations in which 
hyperbolicity is not assumed and we construct asymptotic solutions of arbitrary order. In Sec- 
tion 3, we treat the proposed class of second-order hyperbolic Fuchsian systems, and establish 
the existence of solutions with prescribed asymptotic behavior on the singularity. 

In [7J, we will apply our Fuchsian framework and treat the class of Gowdy spacetimes, and 
this will allow us to recover (and slightly generalize) earlier results by Rendall and collaborators, 
while providing a more direct approach. Although conceptually similar, the proposed analysis 
based on second-order equations lead to a simpler description and provides a definite advantage 
for the applications. Moreover, as we will demonstrate in [7J. the approach introduced in the 
present paper can be cast into a discretization scheme and allows us to numerically and accurately 
compute solutions to the initial value problem. Our theory also turns out to be robust enough to 
extend to matter models, especially to the Einstein-Euler equations [3] [TOJ [13 [H] , as discussed 
in@. 

2 Second-order Fuchsian systems 
2.1 Terminology and objectives 

In this section, we rely mainly on techniques for ordinary differential equations (ODE's). The 
main theory of interest developed in the next section (Section 3) will require an hyperbolicity 
assumption which is not yet made at this stage. The main purpose of this section is to present 
some important terminology and concepts within the simple framework of ODE's, and to point 
out some difficulty arising with singular equations. 

First of all, t > denoting the time variable, the operator 

D := td t 

will often be used, rather than the partial derivative dt- Indeed, the weight t is convenient to 
handle asymptotic expansions near the singularity t — 0. Occasionally, we write D t instead of 
-D, especially when several time variables are involved. 
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Definition 2.1 (Second-order Fuchsian systems). A second-order Fuchsian system is a system 
of partial differential equations of the form 

D 2 u{t,x) + 2A(x) Du(t,x) + B(x) u(t,x) = f[u](t,x) (2.1) 

with unknown function u : (0,5] X U — > R™ (for some S > and interval U C M.), where the 
coefficients A = A(x) and B = B(x) are diagonal, n x n matrix-valued maps defined on U , and 
the source-term f = f[u](t,x) is an n-vector-valued map of the form 

f[u]{t, x) := f{t, x, u, Du, d x u, d x Du, d k x u, d k x Du), (2.2) 

for some integer k > 0. 

We assume that the coefficients A and B do not depend on t and we derive our results 
under this assumption. The generalization to time-dependent coefficients does not bring essential 
difficulties. The assumption that the matrices A, B are diagonal is not a genuine restriction, if 
the system with arbitrary matrices A, B can be recast in diagonal form, i.e. if A, B admit a 
common basis of eigenvectors. Under this condition (which is always satisfied in the applications 
of interest in general relativity), the system is "essentially decoupled" since the coupling takes 
place in terms of non-leading order, only. 

We denote the eigenvalues of A and B by . . . , aS n > and . . . , b^ n \ respectively. When 
it is not necessary to specify the superscripts, we just write a, b to denote any eigenvalues of 
A,B. With this convention, we introduce: 

Ai := a+ \ja? - b, A 2 := a - \j a? - b. (2.3) 

It will turn out that these coefficients describe the expected behavior at t = of general solutions 
to (gZD) . 

In Dcfinition l2.il the assumption that U is a one-dimensional domain makes the presentation 
simpler, but most results below remain valid for arbitrary spatial dimensions. For dchnitcncss 
and without much loss of generality, we assume throughout this paper that all functions under 
consideration are periodic in the spatial variable x and that U is the periodicity domain. All 
data and solutions are extended by periodicity outside the interval U. 

The left-hand side of (|2.ip is referred to as the principal part of the system. The reason for 
incorporating certain lower derivative terms in the principal part is that we expect these terms to 
be of the same leading-order at the singularity t = 0. In contrast, the source-term is anticipated 
as negligible in some sense there, see below. Observe that, at this level of generality, there is some 
freedom in bringing terms from the principal part to the right-hand side, and absorbing them 
into the source-function / (or vice- versa). This freedom has several (interesting) consequences, 
as we will discuss later on: roughly speaking, some normalization will be necessary later, yet at 
this stage, we do not fix the behavior of / at t = 0. 

We are mainly interested in solving a singular initial value problem associated with (12. 1|) . 
with data prescribed on the singularity t = 0, in a sense made precise later on. The fundamental 
question is, of course, to determine conditions on the data and coefficients ensuring existence 
and uniqueness of a solution u. It will turn out that the behavior at t = cannot be prescribed 
arbitrarily, but is tight to the value of the coefficients Ai and A2 defined in (|2.3p . Indeed, we will 
specify the behavior of solutions at t = 0, in terms of freely specifiable functions (the data on the 
singularity), and derive an asymptotic expression of arbitrary order providing the asymptotic 
form of general solutions. 
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2.2 The case of linear ODE's depending on a: as a parameter 
Explicit formula 

We begin our investigation of the second-order Fuchsian systems (|2.ip by treating the case 
/ = w(t, x) for some given function w. We are led to consider a family of scalar ordinary 
differential equations that are completely independent from each other; recall that the matrices 
A and B of the principal part are diagonal. Without loss of generality, we thus assume that 
n = 1 throughout the present section. In turn, the spatial variable x is treated as a parameter 
which we do not need to write explicitly yet. As we will see, it is instructive to express the 
general solution in this elementary case, as we now do. 

Consider the following singular, inhomogeneous, singular ordinary differential equation 



D 2 u(t) + 2aDu(t) + bu(t) = w(t) 



(2.4) 



with unknown u = u(t), where w = w(t) is a given locally integrable function. Further integra- 
bility of w near t = will be imposed shortly below. Recall that a and b are constant in t and, 
Ai, A2 were defined in (|2.3[) . We begin with a formal result and the convergence of the integral 
terms will be discussed rigorously later. 

Proposition 2.2 (Linear second-order Fuchsian ODE. Formal version). General solutions of the 
inhomogeneous singular ordinary differential equation (|2.4p are given by 



u(t) 



u*t a lnt + u** t a + 



w(t/OC 



-a-1 



ln(d(, 



u* t 



-Ai 



f— A2 



(2.5) 



Ai —X'2 



m) (c 



in which u* and m** are prescribed data. Alternatively, one can write 

ri 



u(t) 



u*t a mt 

M* t~ Xl + U** t 



r 



w{s)s° 







In - ds, 

s 



1 



A1-A2 



-A 3 



w(s)s 



A2-1 



ds — t 



-A, 



w(s)s 



Ai-1 



75|) in the form 
a 2 = b, 

a 2 f b. 



Proof. In the rescaled time variable r\ := — hit, equation (|2.4p has constant coefficients, indeed 

u (ri) — 2au {rj) + bu(r]) = w(rj), (2-6) 

where u(rj) := u(e~ ?) ) and w(rj) :— w(e'~ ri ) and the prime ' denotes a derivative with respect to 77. 
This is nothing but a linear harmonic oscillator equation with friction term —2au' and forcing 
term w. The singularity of (|2.4p at t — corresponds to the singularity at infinity 77 = 00. 

First, we seek for general solutions of the homogeneous equation for w = 0. From the ansatz 
u = e Ar ' we obviously get the roots defined in (|2.3p . The general solution of the homogeneous 
equation is thus 



u(r)) 



— u* rj e av + u** e° 



m* e 



Ai?) 



o?=b, 
a 2 ^b, 



where u* and are constants with respect to rj (and the negative sign is chosen for convenience 
in the following discussion). 
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Particular solutions of the general inhomogeneous equation (|2.6I) are easily constructed by 
the Duhamel principle. Let u := u{rj) be the solution of the homogeneous equation for vanishing 
data w(0) = and u'(0) — — 1. For any given w, the function 

r° ~, , ~ 

= / u(t) w (rj — t) dr 

J — oo 

is a particular solution of the inhomogeneous equation. (This is true only formally at this 
stage, since we have not yet checked under which conditions the integral exists and can be 
differentiated.) Hence, the general solution u of (12.61) reads 

!— w* 7/ e ar > + it** e ar i — Te aT w(i] — r)dr, a 2 = b, 
u* e Xir > + u»* e X2r > 

Returning to the original time t > 0, we find 

f lnt + M «t- a -/^(-lnC)C- a w(ln(CA))(-l/C)dC, a 2 = b, 

{ ~xihr 2 IL (C Xl C X2 ) ®(in(C/t))(-i/C)dC, « 2 * 6, 

and by substituting r = — In £, this concludes the proof. □ 
The spatial coordinate x as a parameter 

For the later discussion, it is convenient to write the spatial variable x explicitly as a parameter 
now. Define r(x) := y a(x) 2 — b(x) which might be real or imaginary dependent on the values 
of the coefficients. If there are points xq € U so that r(xo) = and other points x% € U with 
r(a;i) 0, then we will renormalize the coefficients u*(x) and «**(x) in (|2.5p as follows. In order 
to obtain a continuous transition from the non-degenerate case r ^ to the degenerate case 
T = 0, let us first rename the coefficient functions for the case a 2 ^ b in (|2.5p to u* and m**. 
Now if we set 

A , v M*(x) - M**(x)/r(x) „ U«(x) +U^(x)/T(x) 

u*(x) = , u**{x) = , (2.7) 

and choose u*(x), 7i**(x) to be, say, continuous functions, then the function determined by the 
two leading terms in (12. lip is continuous in x for all £ > even at x — xo, provided T is 
continuous. Indeed the full general solution u(t, x) of (|2.4p is continuous in x for all t > in this 
case. 

In view of Proposition 12.21 it is natural to define the solution operator H associated with a 
source function w — w(t, x) by 

t 



t~ a ^ w{s,x)s a(x) - 1 \n-ds, (a(x)) 2 = b(x), 

Jo t s 

?—^(t- x ^ f w{s)s M*)-i ds (2.8) 
Ai(x) - A 2 (x) V J Q 



(H[w])(t,x) 



It represents the solution of (|2.4p for the choice it* = u*, = 0, at least on the formal level so far. 
According to the previous discussion, the non-degenerate case a 2 ^ b in the definition converges 



t~ Al(a;) J w(s,x)s Xl(x) - i y {a{x)f ^b{x). 
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to the degenerate case a? — b at degenerate points continuously, provided the coefficients are 
continuous, and vice versa. 

Fixing some 5 > 0, we now assume that the source w in (|2.4p belongs to C lxm ((0, 5} x U), that 
is, w is /-times continuously differentiable with respect t and m-times continuously differentiable 
with respect to x on (0,<5] x U C K 2 . Here, I and m are non-negative integers. Moreover, we 
assume that the coefficients a and b of the equation are C rn (U). In this case, the general theory of 
ordinary differential equations implies that the solution u(t, x) of (I2.4j) depends as a C m function 
on x. Hence, if tt[w] rigorously represents a particular solution of (|2.4[) . then i) the function 
H[w](t, x) is in C m (U) with respect to x for each t > 0, ii) we can take the spatial derivatives 
under the integral, iii) each spatial derivative of H[w](t,x) converges from the non-degenerate 
to the degenerate case at degenerate points as a continuous function, and vice versa. 

Behavior near the singular time 

Now we go beyond a formal derivation and determine precise conditions on w under which the 
integrals in (I2.5|) make sense and (|2.5p provides actual solutions of (|2.4[) . We use here the notation 
•ft for the real part of a complex number. 

Proposition 2.3 (Pointwise properties of the solution operator H). Fix some S > 0, a compact 
set K C U, and l,m > 0, and let w be a function in C lxm ((0,5] x U), and a, b e C m (U). In 
addition, suppose that w satisfies the following asymptotic conditions: there exists a constant a 
such that 

a > -3?A 2 (x), x e K 

and, for all < p < I, < q < m, 

sup\D p d%w{t,-)\ = 0{t a ). 

K 

Then, the operator H given by (|2.8p is well-defined and, ifl>2, it provides a particular classical 
solution of (|2.4[) . Moreover, we have H[w] G C (xm ((0,<5] x K) with, for all sufficiently small 
e > 0, 

sup \D p d q x H[w]{t,-)\ = 0{t a ~ e ). 

K 

In addition, for 1 < p < I one has 

{D p t H[w]){t, x) = D^wis, x)s a{x) ~ 1 (a(x) In | + l) ds (2.9) 
when a 2 (x) — b(x) and, otherwise, 

(D P H[w])(t, x) = ( - \ 2 {x)t-^ f D^wis, x)s x ^~ 1 ds 

v f2.im 
rt 

+ Ai(x)^ Al ^ / D p s - 1 w{s,x)s x ^ x) - 1 ds 
Jo 

We note that we are allowed to choose e — in the previous proposition only if a(x) and b(x) 
are constants in space or for q = 0. The constant e > is necessary in order to control logarithms 
which arise when spatial derivatives are taken of functions involving spatially dependent powers 
oft. 
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Proof. The proof is quite direct and we only show the derivation of the formulas (|2.9p and (I2.10j) . 
For definiteness, we treat the case a 2 (a;) 7^ b(x) only, since the other case is treated similarly 
and the transition between the two cases can be obtained by appropriate limiting procedures. 
Consider the expression for H in (12.5[) . Under our regularity assumptions, all derivatives are 
calculated by differentiation under the integral sign, and we obtain 

1 f°° 

(ITH[w])(t,x) = ^—^ J D*w{t/Z,x) (r^" 1 -C Xl - l )dd. 

As done earlier, we introduce the new variable s := t/C and observe that D t = tdt = sd s =: D s 
for fixed £. Hence, we obtain 

(D*H [w])(t, x) = [ BP w (a, x) (C^ 1 - C Al_1 ) -ds 

M - A 2 Jo s 



= - — / d s {D p s - 1 w){s,x) (r x *s Xa -r Xl s Xl )ds 

Ai — A 2 J 

= x 1 _ x ^' A2 (^D p s - 1 w(s,x)s X2 \ t - X 2 D p ~ 1 w{s,x)s X2 - 1 d^j 

-t- Xl (^D p s - 1 w(s,x)s Xl \ t Q - \i J D p s ~ 1 w(s,x)s Xl - 1 ds 

where we used integration by parts. Now, in view of our regularity assumptions, we conclude 
that all terms here have a limit when t — > 0. All terms, except for the main integrals, either 
vanish or cancel. □ 



2.3 Asymptotic solutions of second-order Fuchsian systems 
Asymptotic data and canonical expansions 

We now return to the non- linear problem (|2.ip . We are going to identify the "canonical" asymp- 
totic behavior (at t = 0) of general solutions to second-order Fuchsian equations using first 
heuristic arguments, and we determine a "canonical" expansion. Such an expansion involves cer- 
tain free functions, interpreted as "data on the singularity" and allows us to formulate a singular 
initial value problem. 

The basic understanding of the behavior of Gowdy solutions at t = is obtained from the 
BKL conjecture, and the idea (which we explain in further detail in the second paper [7]) is 
to neglect spatial derivatives in the evolution equations while solving the remaining ordinary 
differential equations at each spatial point x. This approach allows to identify the leading order 
terms in the expansion of the solution at t = 0. Indeed, the Gowdy equations turn out to be 
of second-order Fuchsian form, as shown in the second paper. This suggests that we use similar 
arguments for the derivation of solutions of general second-order Fuchsian equations. 

According to the above heuristics, the behavior of general solutions to (|2.1j) should be driven 
by the principal part of the PDE's with coefficients evaluated at t = if the source-term satisfies 
certain "decay properties" at t = 0, as we will discuss later. More precisely, the two leading 
terms in the expansion of solutions at t = should be determined by the homogeneous equa- 
tions obtained by setting the right-hand side / in (|2.1j) to zero. According to the discussion in 
Section 12.21 any solution u to the second-order Fuchsian equations (|2.ip should hence have an 
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expansion of the form 

*£° (») t~ Q<l) (x) In i + nil (x) t- a{%) (*) + 0(^ q(i) ) , 

I ( i )( a; ))2 =6(0(3 



(i,z) 



< 



,00/ 



CO, 



-\<£\x) 



(2-11) 



for each i = 1, . . . , n and each x £ U . Here, the functions and u** are prescribed, and offl > 
are real constants. The meaning of the Landau symbols O in this context will be made precise 
later; at this stage of the discussion they have to be understood "intuitively" as representing 
terms of higher order in t at t = 0. In the case of a transition from the non-degenerate case to 
the degenerate one or vice versa, the renormalization given by (|2.7[) is necessary and will always 
be assumed. 

At this stage of the discussion, we clearly see the dependence of the expected leading-order 
behavior at t = on the coefficients of the principal part of the equation. If the roots Ai and 
A2 are real and distinct, i.e. if a 2 > 6, we expect a power-Jaw behavior. In the degenerate case 
Ai = A2, i.e. if a 2 — b, we expect a logarithmic behavior. Finally, when Ai and A2 are complex 
for a 2 < b, the solution is expected to have an oscillatory behavior at t = of the form 

u(t, x) = t~ a ^{u* cos(A/(x) log*) + u** sin(A/(a;) log*)) + . . . 

for some real coefficient functions u*[x) and u**(x); note that in this case, Ai = A2 = a + i\j 
with Xj := \/b 2 — a. 



Relevant function spaces 

Consider any second-order Fuchsian system of the form described in Definition 12. 1[ with coef- 
ficients a, b, Ai,A2 satisfying (|2.3[) . To simplify the presentation, we restrict attention to scalar 
equations (n — 1) and shortly comment on the general case in the course of the discussion. 

Fix some integers I, m > and constants a, S > 0. For w 6 C ((0, 5], H m (U)), we define the 
norm 

i \ 1/2 

t mX2ix) ~ a) \d q x D p w(t,x)\ 2 dx\ , (2.12) 

and denote by Xs,a,i,m the space of all functions with finite norm ||io||<s )a j )m < 00. Throughout, 
H m (U) denotes the standard Sobolev space and we recall that all functions are periodic in the 
variable x with U being a periodicity domain. To cover a system of n > 1 second-order Fuchsian 
equations, the norm above is defined by summing over all vector components with different 
exponents used for different components; recall that each equation in the system will have a 
different root function A2 and we allow that a = (a^ , . . . , a*")) is a vector of different positive 
constants for each equation. The constant 5, however, is assumed to be common for all equations 
in the system. With this modification, all results in the present section remain valid for systems 
of equations. 

Throughout it is assumed that 3?A2 is continuous and it is then easy to check that {Xs, a ,l,mi II 1 
IU,a,z,m) is a Banach space and that the following property holds. 

Lemma 2.4 (Approximation by smooth functions). Given any w € Xs,a,l,m> a-nd constant e > 0, 
there exists a sequence (w n ) € Xs >a ,l,m H C°°((0, S] x U) such that 

lim \\w - W n \\s, a - e ,l,m = 0. (2.13) 
n—>oo 



\W\\s,a,L: 



sup 

0<t<6 
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The functions w n above are taken to be periodic in space, for instance: 

/oo poo -I 
/ t x ^w(s,y)kjln-)k v (x-y)-dsdy. (2.14) 
oo JO \ 1/ S 

Here, k v : K — » M+ is any smooth kernel supported in [—77,77], satisfying j R k v (x)dx = 1 for 
all positive 77. We note that the constant e > in (|2 . 13[) is introduced in order to guarantee 
uniform convergence on (0, S] in the case where w has no limit at t = 0, being just bounded and 
continuous. 

In terms of the spaces -X^a^m, we define in a mathematically precise way a notion of canonical 
expansions and asymptotic data as follows. 

Definition 2.5. Consider a second-order Fuchsian equation (|2.1j) with continuous coefficients 
a, b, Ai, A2. Suppose that v and w are functions related as follows: 

, {u*(x)t- a ( x hnt + u**(x)t- a W+w(t,x), (a(x)) 2 =b(x), 

{,X) \u*(x)t- x ^+u**(x)t- X2 W+w{t,x), (a(x)) 2 ^b{x), { ' ' 

for some prescribed data and it** e H m (U), where m! is some non-negative integer. Then, 
one says that v satisfies a canonical two-term expansion with asymptotic data it* and it** 
and remainder w, provided w g Xs. a ,i, m for some constants 6, a > and non-negative integers 
I, m. 

If the coefficients of the equations are such that there is a continuous transition between the 
two cases in (|2.15p . then the asymptotic data functions it, and are renormalized by (|2.7p . 

2.4 An existence result for second-order Fuchsian ODE systems 

An important property of the ODE solution operator H defined in (|2.8[) is derived now. 

Proposition 2.6 (Continuity of the ODE solution operator H). Pick up any constants S > 0, 
a > and any integers I > 1, m > 0. Then, for each e > 0, the operator H defined in (|2.8I) 
extends uniquely to a continuous linear map Xs, a ,l—l,m ~^ Xs. a ~e,i.m, arid there exists a constant 
C e > (independent of 6 provided S is sufficiently small), so that 

\\H[w]\\ s < C e 6 e \\w\\ s t a,l—l,m 1 

(2.16) 

for all w € Xs, a ,l-l,m- 

We stress at this stage that for all Z > 2 the extended solution operator H indeed provides 
the general solution to the Fuchsian equation: given any function g € Xs^ a ,i-i,m, the function 
w :— H[g] £ Xs. a -e,i,m satisfies the second-order equation 

D 2 w + 2a Dw + bw = g, 

as equality between functions in the space Xs >a -n,i-2.m- 

Proof of Proposition \2.b\ We restrict attention to the case a 2 7^ b, the proof of the other case 
is completely similar and the transition between the two cases can be understood as a limiting 
process. Consider first a function w € Xs^ a ,i-i,m H C°°((0,<5] x U). We have 

\\H[w]\\l a _^ m = sup / t 2 ^-^Y.H\ d " DPH ^\ 2d ^ 
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which is finite thanks to Proposition 12.31 for all e > 0. The term p = is expressed explicitly by 
means of (|2.8p . while the case p > is treated in (|2.9p and (|2.10j) . Suppose (i, x) e (0,(5] x £/. 
For convenience, we introduce functions w\ and w 2 so that 



For p — 0, these functions are 

Wi (i, a;, s) : 
w 2 (t, x, s) : 



d q (w 2 (t, x, s) — w\ (t, x, s))s 1 ds. 



1 



Ai(x) - \ 2 (x) 
1 

Ai(a;) - \ 2 {x) 



w(s,x)t~ Xl{x) s Xl( - x \ 
w{s,x)t- X2{x \s x ^ x \ 



while for p > we only need to substitute w by D p 1 w and complement the expression by factors 
whose particular form is not relevant for the following. With this, we get 



t 2 
d q (w 2 (t, x, s) — wi(t, x, s))s~ 1 ds 



\d q x D p H[w](x)\ 2 



The first term (and in the same way the second one) is handled via Cauchy-Schwarz's inequality 
(for some constant r\ > 0) 







i 




< 2 


/ d%w 2 s 1 ds 


+ 2 


/ d%w\s x ds 




Jo 




Jo 



f |a> 2 s"Vs = / \d q x w 2 \s-^ +r >^ 2 s- {1 ^' 2 ds 
Jo Jo 



< 



ft 



2 a- 1 -" da 



1/2 



1/2 



-f f mwrfs-^da 
V Jo 



1/2 



Hence, we have 

t 2 ^- a+ ^\d q x D p H[w]{x)\ 2 dx 



< 



(/ L tmX2 ~ a+£)tV ( d >i) 2s ~ 1 ~ Vdsdx + 1 1 t^-^h^diwrfs-^dsdx). 



V ^Ju Jo 



(2.17) 



The first term on the right side of this inequality can be written as 

t m>>2-*+e) t V( d ci W2 f s -l-V dsdx 

£ (J ( s 2 ^ x *- a Hd q w 2 t x \s- x *) 2 s 2 ^h 2 ^- e 



u Jo 



g \ 2(a-ei)+2iSA 2 



< 



'■') G) 

-t 2 ^ sup ( [ ( S 2 ^-- a \d q w 2 t x * S - x *) 2 S 2 ^h 2 ^)dx). 
V 0<s<t \Ju v 1 J 



dx)t 1 i S - 1+r ids (2.18) 



11 



Here we have assumed that the constant e\ satisfies e% < a. The significance of the terms s 2 ^ - '') 
and t 2 ( e_Cl ) becomes clear in a moment. For the second term on the right side of (I2.17p . we get 
similarly 

t 2 (^2-a+e) t r,^ d q w ^2 s -l-n dsdx 



(7 JO 



< V> sup f / U^- a \dlw^s-^fs^-^-^\dx 

V 0<s<t \Ju v 



(2.19) 



where in the last step we used that SRAi > SRA2. Now, for p = 0we compute, 

Ai — A2 

etc.; analogous formulas hold for W2 and for p > 0. So the terms d^w\t Xl s~ Xl incorporate spatial 
derivatives of w of all order lower or equal to q, and all orders lower than q are multiplied with 
logarithmic terms in t and s. Hence, in order to guarantee that the suprema in (|2.18p and (|2.19[) 
are finite for each t > 0, we must choose e\ > r\. Moreover, the suprema are uniformly bounded 
for all t € (0, 8], if e > e\. With these choices, we can rearrange all terms, introduce a finite 
constant C, which is independent of 8 if 8 is small, as in the hypothesis, and hence obtain 

\\H[w]\\ S , a -e,l, m < C8 n \\w\\s, a ,i~i, m , 

for all w € Xs t a,l—i,m n C°°((0,8] x U). Since this inequality holds for all 77 smaller than e, it 
also holds in the limit 77 — > e. The constant C can hence be adapted so that (|2.16[) follows for all 

»a W -i !m nr((o,i] xu). 

Now let w be a general element in Xs^ a ,i-i,m and choose a positive eo with eo < £■ We set 
e := e - e . According to Lemma H3 there exists a sequence (iu n ) C Xj, nC°°((0,<5] x 17), 

so that limn-j.oo ||io — u; n ||,5 )Q ._ e0j ;_i )m = 0. Our results for the smooth case show that (H[w n ]) 
is a Cauchy sequence in Xs. a -e -7,i-i, m = Xs. a -e,i-i, m , and we can denote its limit element by 
H[w] € Xs^ a - et i-i t m- In this way, we define the extension of H to the whole space Xg iCt i—i tfn . It 
is then straightforward to see that the limit of the estimate (|2.16[) leads to the claimed estimate 
for the full space Xs^ a -e Q ,i-i,m- d 

Consider a Fuchsian equation (|2.1[) with right-hand side / of the form (|2.2j) . Let v be a 
function in the form of Definition 12 . 51 with remainder w € Xs, a ,i,m and prescribed data u*, u** G 
i7 m ([/) for some 5 > 0, a > and i, m, m' > 0. Suppose that m, m' are sufficiently large. If 
necessary, the derivatives in (|2.2p are understood in the sense of distributions and we set 

F[w}(t,x) := f[v}(t,x) (2.20) 

and finally 

G:=HoF. (2.21) 
We are now in a position to define an iteration sequence based on this operator G. 
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Proposition 2.7 (Iteration sequence). With the same notation and assumptions as in Propo- 
sition \2.b\ let k be the number of spatial derivatives in f according to (|2.2[) and mo another 
non-negative integer. Suppose that, for given asymptotic data, the operator F satisfies the fol- 
lowing regularity assumption (for some 6q > 0) 

F : Xs, a ,l,m —> -^5,a+e ,;-l,TO-fej (2.22) 

for an integer I > 1, and all non-negative integers to with k < m < tuq. Let G be the operator de- 
fined in (|2.2ip . Then, given any w\ € Xs^ a ,^ mo , the (in general finite) sequence (wj) determined 
by 

Wj+i = G[uij\, for all integers j £ [l,mo/fc] (2.23) 
is well-defined and, moreover, uij+i € Xs. a ,i.m -jk- 

Consider a sequence (vjj ) defined in the above lemma. It determines a sequence (vj ) for fixed 
asymptotic data according to (|2.15[) . and all Vj satisfy the canonical two- term expansion. For 
hi ^ 2, the function uij € Xs. a ,i. m satisfies the second-order equation 

D 2 Wj + 2aDvjj + buij — F[wj-i] 

as equality in the space Xs, a ,i—2,m- The sequence (uij) has infinitely many elements if mo — oo ot 
if k = 0. In the latter case the second-order Fuchsian equation is a system of ordinary differential 
equations (with the spatial variable parameter). In both function w is a fixed 

point of the iteration sequence if and only if the associated function v is a solution of (|2.1j) . A 
fixed point theorem for the ODE case is as follows. 

Theorem 2.8 (Existence of solutions to second-order Fuchsian ODEs). Under the assumptions 
as in Proposition \2.1\ with k — 0, suppose additionally that for given asymptotic data, the operator 
F satisfies the following Lipschitz continuity property; for each r > and eo > arising in 
(|2.22l) . there exists C > independent of 8, so that 

\\F[w] - F[w\\\ 5ta+€Q ,i-x,m < C\\w- w\\ s , a ,i, m (2.24) 

for all w,w G B r (0) C Xs ia ,i,m- Then, given any initial data w\ £ Xs t0lt i iTn and provided S > 
is sufficiently small, the iteration sequence (|2.23p converges to a unique fixed point w € Xs^ a .i,m- 

Proof. Our previous results imply 

\\G[w] - G[w\\\s, a ,l, m < C5 n \\w - w\\s,a,l,m 

for a constant C > 0, provided w,w £ Bk{0) C Xs,a,l,m- Hence, for sufficiently small 8, the 
operator G becomes a contraction. The convergence of the iteration sequence follows if we can 
guarantee that Wj € Br(0) for a sufficiently large K. This, however, is the case since the sequence 
(uij) is a Cauchy sequence thanks to the contraction property of G and, hence, is bounded. □ 

Condition (|2.24l) guarantees convergence of the sequence (wj) in the ODE case. The approach 
of this section is not sufficient to cover the PDE cases of interest and, in general, this condition 
does not hold if k > 0. Moreover, the iteration sequence constructed above is only finite if too 
is finite. We can expect that in typical applications, too is infinite if the asymptotic data are 
smooth and, say, w± = 0. Still, this does not lead to an existence result except for the analytic 
case, see [14]. Well-posedness for second-order Fuchsian PDE's in a larger than the analytic 
class will be addressed in Section [3] after we make the additional assumption that the Fuchsian 
equations are hyperbolic. 
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2.5 Asymptotic solutions of arbitrary order 

The iterative sequence (wj) has useful asymptotic properties. In order to simplify the discussion 
in this section, we assume, instead of (I2.2p . that / has the form 

f[u](t,x) := f{t,x;u,Du,d x u,d x Du,d 2 x u) (2.25) 

and that it is a polynomial in all of the arguments involving u with coefficients which are smooth 
and spatially periodic on (0, S] x U. The operator F is defined, in the same way as was done 
earlier, from given asymptotic data u* and it*». We henceforth assume in the following that 
u*, it*» £ H mi (U) for some non- negative integer mi and that (I2.22p holds for k = 2 and mo = mi. 

Definition 2.9. A function v satisfying the canonical two-term expansion with given asymptotic 
data is called an asymptotic solution of order 7 > to the system (|2.ip provided the residual 

R[w] := L[w] - F[w] 

belongs to ^C5, 7 +j}A 2 ,o.O; * n which the following notation 

L[w]:= D 2 w + 2aDw + bw, (2.26) 

is used for the principal part of (|2.1[) and one assumes that w £ Xs, a ,i, m with l,m>2. 

In this definition, we use the obvious generalization of the spaces Xs,si,l,m to spatially depen- 
dent exponents a. Note that L[w] — L[v] if W and v are related as in (|2 . 1 5[) . When I = 1 or 
m = 1, one needs to reformulate the operator L (hence R) in a weak form, as we will explain 
below. 

Theorem 2.10. Suppose that the operator F satisfies the conditions stated earlier for given 
asymptotic data, and consider the iteration sequence Wj £ Xg a i mi _ 2 (j-i) for j > 2 given by 
(|2.23[) with u>i = 0. Then, for any constant k < 1 the sequence has the property 

Wj+l — Wj £ ^(5,a+(j-l)Ke ,;,mi-2j 

for 1 < j < mi/2. Moreover, the residual satisfies 

R[wj] £ ^5,a!+(j-l)Ke ,i-l,mi-2(j-2)j 

and hence, Wj is an asymptotic solution of order 7 = — SRA2 + ot + (J — 1)k£o f or 2 < j ■ 

This establishes that the order of the asymptotic solution Wj is an increasing function in j. 
Hence, the functions Wj can be interpreted as approximations of actual solutions of increasing 
accuracy at t = 0. 

Proof. In order to show the first relation, we proceed inductively and start with j = 1. We need 
to show that w 2 £ ^i, a ,i.mi-2 which is true by Proposition 12.71 Next, we suppose that 

w j ~ w j-l G ^<5,a+(j'-2)«eo,(,mi-2(j-l)i 

has already been shown for a given integer j £ [2,mi/2 + 1]. We take the difference of the 
equations for u>j+i and Wj, and the linearity of the operator H implies 

w j+1 - wj = H[F[w,} - F[wj-i]]. (2.27) 
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Since / of the form (|2.25p depends smoothly on its arguments by assumption, the mean value 
theorem implies the existence of a matrix- valued function M of the form 

M[uij, Wj-i](t, x) := M(t, x; Wj, Dwj, d x Wj,d x Dwj,d x Wj, 

Wj-i,Dwj-i,d x Wj-i,d x Dwj-i,d^.Wj-i) 

depending as a polynomial on all arguments involving Wj and with the property that first 

M[vjj,Wj-i] G Xs t snx 2 -a+eo,l-l,m 1 -2{j-2)> 

and that second 

F[v)j] - F[wj-i] = M[wj,Wj-i] ■ AV[wj,Wj-i}. (2.28) 
Here, the vector-valued function AV is defined as 

AV[wj,Wj-i] := (wj -Wj-i,Dwj - Dwj-i, d x Wj — d x Wj-i, 
d x Dvjj - d x Dwj-i, . . . ,dl<Wj - d x l Wj- 1 ) T . 

We have, 

AV[Wj, Wj-l] e ^5,Q+(j-2)K£o,/-l,mi-2(i-2)- 

All this implies that 

F l w j] - F[lVj-l] G -X" 5>Q , + ( J -_2) reeo+eo ,i-l, mi -2(j-2)- 

Then, Proposition 12.61 yields 

Wj+l — Wj G ^5,a+(j-l)Ke ,i,mi-ib(j-2)- 

At this point we see the significance of the requirement k < 1. Namely, we must choose the 
constant e in Proposition 12.61 as e = eo(l — «). 

The second claim of the theorem is now an immediate consequence. We write the system, 
which determines Wj , in the form 

- F[wj] = F[wj-i] ~ F[wj]. 

Again, we can write the right side as above. We conclude from the previous results that the right 
side is in Ar^ ct+ ( J _ 2 ) Keo + £o ./_i. mi _2(j-2) f° r 2 < j. Since the left side equals R[wj], the result 
follows. □ 

3 Second-order hyperbolic Fuchsian systems 
3.1 Assumptions and basic definitions 

From now on we focus on hyperbolic second-order Fuchsian equations in the sense of Definition l3.il 
below — which form a special case of general second-order Fuchsian equations. Our aim of this 
section is to establish a well-posedness theory for the (singular) initial value problem when data 
are prescribed on the singularity. 

Definition 3.1 (Second-order hyperbolic Fuchsian systems). A second-order hyperbolic Fuchsian 
system is a set of partial differential equations of the form 

D 2 v + 2ADv + Bv-t 2 K 2 dlv = f[v], (3.1) 
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in which the function v : (0, 8] x U — > WL n is the main unknown ( defined for some 8 > and some 
interval U), while the coefficients A = A{x), B = B{x), K — K(t,x) are diagonal nx n matrix- 
valued maps and are smooth in x € U and t in the half-open interval (0, 8], and f — f[v](t, x) is 
an n-vector-valued map of the following form 

f[v](t,x) := f(t,x,v(t,x),Dv(t,x),tK(t,x)d x v(t,xfj. (3.2) 

As earlier, we assume that all functions are periodic with respect to x and that U is a period- 
icity domain. Further restrictions on the coefficients and on the right-hand side will be imposed 
and discussed in the course of our investigation. Hence, hyperbolic second-order Fuchsian sys- 
tems are second-order Fuchsian systems with a particular structure of their right-hand side: in 
our new notation, we have separated the second-order spatial derivatives from other terms in the 
right-hand side / and incorporate them into the principal part of the equation, which now reads 

L:=D 2 + 2AD + B -t 2 K 2 d 2 x . (3.3) 



This is a linear wave operator for t > and, indeed, (I3.ip is hyperbolic provided K satisfies 
(positivity) conditions given below. Later on we will construct solutions where the first three 
terms of the principal part are of the same order at t — and dominant as in the previous section, 
while the second spatial derivative term is assumed to be of higher order at t = 0. Hence, we 
expect the same phenomenology at t = as earlier, and the only significance of the new term in 
the principal part is that it allows to derive energy estimates. 

The eigenvalues of the matrix K are denoted by and, in the scalar case (or when there is 
no need to specify the index), we simply write k. These quantities are interpreted as characteristic 
speeds. Throughout this section, we assume that they have the form 

with fi {l) : [/-> (-1 , oo), v^' : [0, 8] x U — > (0, oo) smooth functions. 

In particular, we assume that each derivative of has a finite limit at t = for each x € U. 
The motivation for these assumptions will become clear in the forthcoming discussion. Note we 
allow for the characteristic speeds to diverge at t = 0. At a first glance, this appears to conflict 
with the standard finite domain of dependence property of hyperbolic equations. Recall that 
for the standard initial value problem of hyperbolic equations, the solution at a given point is 
determined by the restriction of the data to a bounded domain of the initial hypersurface; this 
is a consequence of the finitcness of the characteristic speeds. A closer look at the requirement 
fi{x) > —1, however, indicates that the characteristic curves are integrable at t = and, hence 
that the Unite domain of dependence property is preserved under our assumptions. 

In order to simplify the presentation and without genuine loss of generality, we restrict the 
discussion now to the scalar case n = 1. Consider any second-order Fuchsian hyperbolic equation 
and for each non-negative integer I and real numbers 8, a > 0, define the space Xs^ a ,i by 

l 

p=0 



and introduce the norm 



( 1 2 

U,a,l ■= I ||/|lla, p ,i- p j i / £ ^,o,I- 

\p=0 ) 
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Recall that the spaces Xs^ a ,i,m and norms || • ||$ )a j,7n have been introduced in the previous 
section. It is straightforward to see that the spaces (Xg tai i, \\ ■ \\s, a ,l) have similar properties 
as the previously defined ones. As we will see in the following discussion, it is not possible to 
control solutions of our equations in the spaces Xs, a j directly. It turns out that we must use 
spaces (Xs, a ,i, \\ ■ W'ga i) instead which are defined as earlier, but in the norm ||/|| £ a ; of some 
function /, the highest spatial derivative term d l x f is weighted with the factor t@ +1 . Here f3 is the 
characteristic speed of the equation given by (13.4[) . It is easy to see under the earlier conditions 
that (Xs ta ,i-, || • ll^a i) are Banach spaces. Moreover, for any w € Xs, a .ii the mollified function 
defined by (|2.14l) is an element of Xs. a - e ,i H C°°((0, 8] x U) for every e > 0. Furthermore, the 
sequence of mollified functions in the limit r\ — > converges to w in the norm || • ||^ Q _ £ We 
also note that Xg_ a ,i C Xs_ a ,i- For our later discussion, let us define 

oo 

Xg,a, oo '■= X$ tCt j, 
1=0 

and note that X s>oti00 = f|^ Xg l0t: i. 

For w £ Xs. a ,i, the operator L in (|3.3p is defined in the sense of distributions, only, via the 
following weak form: 

:= f (*)-«/ -Dw(t,x)D(t)(t,x) + (2A(x)-$l\2(x)+a-l)Dw(t,x)<j)(t,x) 

JO JR ^ 

+ B(x) w(t, x)(f)(t, x) + tK(t, x)d x w(t, x)tK(t, x)d x <j)(t, x) 

+ (2td x K(t,x) +d x ^.X 2 (x)K(t,x)t\iit)tK(t,x)d x w(t,x)(f>(t,x)j dxdt, 

(3.5) 

where 4> is an y test function, i.e. a real- valued C°°-function on (0,5] x M. together with some 
T G (0, S) and a compact set K £ K so that (f>(t, x) = for all t > T and x g" K, and each 
derivative of <p has a finite (not necessarily vanishing) limit at t = for every x € U. The 
formula (|3.5[) is obtained as follows. First we assume that w is a smooth function for t > in 
Xs^a.i- We compute L[w] and multiply (|3.3p with t^* 2 ~ a . Then we integrate this expression in 
x on U and in t on [e, 8] for some e > 0, and integrate by parts. The resulting expression, which 
resembles (|3.5p plus a boundary term at t = e, is meaningful for general w € Xs, a ,i and only 
the limit e — > remains to be checked. It turns out that the assumption w G Xs^ a ,i is sufficient 
to guarantee that this limit is finite and in particular that the boundary term at t = e resulting 
from integration by parts goes to zero in the limit. For our later discussion, we note that for any 
given test function <j), the linear functional (£[•],</>) : X$ ta \ — > R is continuous with respect to 
the norm || • ||^ Q 1 . This is the main reason to include the factor t^2(x)-a j n ^ e definition of C. 
Consider now functions u, v, w on (0, 8] x U related by 

v(t,x) = u(t,x) + w(t,x). (3-6) 

In the following, v will stand for a solution of a Fuchsian system. The function u will be called 
the leading-order part and w the remainder of the solution at the singularity at t = 0. In 
agreement with the discussion of the canonical two-term expansion in Definition [531 we will look 
for remainders w in spaces Xs^ a ,i for some a > 0. However, at this stage of the discussion we 
will not yet fix the particular form of the function u and its dependence on the asymptotic data. 
Indeed let us assume that u is some given function. In analogy to our earlier discussion, we 



17 



introduce the operator F as 

F[w](t,x) := f[u + w](t,x). 

If the operator F is a map Xs. a ,i — > Xs_ a fl,w <— > F[w], which we shall assume later on, it is 
meaningful to define its weak form by (for all test functions </>) 

(F[w],(f>) := f f t iRX2( - x ^ a F[w]{t 7 x)(l)(t,x)dxdt. 
Jo Js. 

Definition 3.2 (Weak solutions of second-order hyperbolic Fuchsian systems). Let u be a given 
function and 8, a > be constants. Then, one says that w € Xg l0ti i is a weak solution to the 
second-order hyperbolic Fuchsian equation (|3.8p . provided 

V[w] := C[w] + C[u] - T[w] = 0. (3.7) 

For our later discussion, we note that for any given test function <f> and under our earlier 
assumptions, the linear functional ("P [•],</>) on Xs, a ,i is continuous with respect to the norm 
II ' \\7u i- F° r completeness, we also state here that the classical form of the equation for w 

L[w] = F[w] - L[u], (3.8) 

if v is a classical solution of the equation (|3.1[) . Clearly, if w and u are sufficiently smooth, then 
w is a weak solution if and only if w is a classical solution of (|3.8[) , or equivalently if v given by 
(|3.6[) is a classical solution to the original second-order hyperbolic Fuchsian equation (|3.1[) . 
We arrive at the following important notion. 

Definition 3.3 (Singular initial value problem (SIVP)). Consider a second-order hyperbolic 
Fuchsian equation (|3.8[) with coefficients (a,b,Xi,\2) and characteristic speeds k satisfying all 
the conditions stated earlier. Moreover, choose a leading- order part u. Then, a function v : 
(0,(5] x U — > R is called a solution of the singular initial value problem provided w := v — u 
belongs to Xs^a.i for some a > 0, and is a weak solution to the second-order hyperbolic Fuchsian 
system (|3.7p . 

In particular, we will be interested in the case when u is parametrized by asymptotic data 
in analogy to the canonical two-term expansion. At this stage of the discussion, however, the 
particular form of u is not fixed yet. 

3.2 Linear theory in the space Xs, a ,i- Main statement 

In this subsection and the following one, we study a particularly fundamental case described by 
the two conditions: 

1. Vanishing leading-order part: u = 0. 

2. Linear source-term: 

F[w](t, x) = f (t, x) + h(t, x)w + f 2 (t, x)Dw + / 3 (f, x)tkd x w, (3.9) 

with given functions fo, f\, fi, fo, so that fx, f 2 , fs are smooth spatially periodic on 
(0, 5] x U, and near t = 

sup/„(i,aj) = 0(^), a = 1,2, 3, (3.10) 

x£U 

for some constant /i > 0. 
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We have not made any assumptions for the function fo yet, since in the following discussion this 
function will play a different role than fx, f 2 , fa. Moreover, no loss of generality is implied by 
the condition u = 0, since the general case can be recovered by absorbing L[u] into the function 
/o. 

Under these assumptions, we pose the question whether there exists a unique weak solution 
w of the given second-order hyperbolic equation in Xs t a,i f° r some 8, a > 0. 

Proposition 3.4 (Existence of solutions of the linear singular initial value problem in Xs, a ,x)- 
Under the assumptions made so far, there exists a unique solution w € Xg. a ,i of the singular 
initial value problem for given 8, a > provided: 

1. The matrix 

/»(Ai-A a ) + a ((^Ai) 2 /r? - r,)/2 

N := (mi) 2 /v ~ ri)/2 a td x k - d x $t(Xx - 

\ td x k-d x $l(\i -A 2 )(tfelnt) 5R(Ai-A 2 ) + a 

is positive semidefinite at each (i, x) € (0, 8) x U for a constant 7] > 0. 

2. The source-term function /o is in Xs^ a +e,o for some e > 0. 
Then, the solution operator 

H : X S ,a+e,0 -> Xs,a,l, fo ^ W, 

is continuous and there exists a finite constant C e > so that 

||H[/ ]||?; a)1 <5 e C e ||/o|| M+ e,o, (3.12) 
for all fo. The constant C e can depend on 8, but is bounded for all small 8. 

For reasons that will become clear later on, we call N the energy dissipation matrix. We 

have assumed that a is a positive constant. If, however, a is a positive spatially periodic function 
in C 1 (J7), the definition of the spaces Xs^ a ,k and Xs iCc .k remains the same, and only the (2, 3)- and 
(3, 2)-components of the energy dissipation matrix N change to td x k — d x QR(Xi — A2) +a)(tk hit). 
In the following, we continue to assume that a is a constant in order to keep the presentation as 
simple as possible, but we stress that all following results hold (with this slight change of N) if 
a is a function, and hence no new difficulty arise. 

3.3 Linear theory in the space Xs :a ,i- The proof 

The main idea for the proof is to approximate a solution of the singular initial value problem by 
a sequence of solutions of regular initial value problems. 

Definition 3.5 (Regular initial value problem (RIVP)). Fix to £ [0,8] and some smooth periodic 
functions g,h : U — > R, and suppose that the right-hand side is of the form Q3.9p with given smooth 
spatially periodic functions fo, fx, f%> f?> on [to, 8] x U. Then, w : [to, 8] x U — > K is called a 
solution of the regular initial value problem associated with the regular data g, h if (|3.8p 

holds everywhere on (to, 8] x U and, moreover, 

w(t ,x) = g(x), d t w(t ,x) = h{x). 



X 2 )(tk hit) 
- 1 - Dk/k / 

(3.11) 
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For the regular initial value problem, we indeed assume that /o is smooth just as f\, f% and 
fa. By the general theory of linear hyperbolic equations, the regular initial value problem is 
well-posed, in the sense that there exists a unique smooth solution w defined on [to, 6] for any 
choice of smooth initial data. 

In order to simplify the presentation, we restrict to the scalar case n = 1 for this whole section; 
the general case can be obtained with the same ideas. Choose S, a > and let w £ C 1 ((0, <5] x U) 
be a spatially periodic function. Then, we define its energy at the time i <E (0, 8] by 

E[w](t) :=e- Rt ~' f t 2 ^ x) - a) e[w](t,x)dx, 

Ju (3.13) 
e[w](t,x) :=- ((r,w{t,x)) 2 + {Dw(t,x)) 2 + (tk(t,x)d x w(t,x)) 2 ) , 

for some constants k > 0, 7 > and 77 > 0. For convenience, we also introduce the following 
notation. For any scalar-valued function w, we define the vector-valued function 

w(t, x) := t m2< - x) - a (r)w(t, x),Dw(t, x), tk(t, x)d x w(t, x)), (3.14) 

involving the same constants as in the energy. Then, we can write 

E[w]{t) = l -e-^\\w(t,-)\\ 2 L2(u) , (3.15) 

the norm here being the Euclidean L 2 -norm for vector-valued functions in x. It is important 
to realize that, provided n > 0, the expression sup 0<t <5 \\w(t, ^W^tu) f° r functions of the form 
(|3.14l) yields a norm which is equivalent to || • ||^ Q v thanks to (I3.4p . Therefore, the energy (|3.13l) 

is of relevance for the discussion of functions in Xs t a,i- 

Lemma 3.6 (Fundamental energy estimate for the regular initial value problem). Suppose that 
the source-term is of the form (|3.9p with the conditions (|3.10[) and that the energy dissipation 
matrix p. Ill) is positive semidefinite on (0,(5] x U for given constants a, n > 0. Then, if S > 
is sufficiently small, there exist constants C, re, 7 > 0, independent of the choice of to £ (0,(5], so 
that for all solutions w of the regular initial value problem with smooth regular data at t = to, 
we have 

\\w(t,-)\\ L 2 (U) 

<Ce^ ty -W(\\w(t Q ,-)\\ L2(u) + ^Vi^-°</o0v)l|£»(iO<fe)> 
for all t G [to , <5] . 
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The role of the matrix N in (|3.11[) in the proof of this result motivates the name "energy 
dissipation matrix" . Moreover, this results demonstrates the importance of the assumption 
f3(x) > —1 in (I3.4[) . Namely, if (3(x) < — 1 at a point x g U, then for any choice of a and 77, 
the matrix N would not be positive semidefinite for small t at x. While the energy estimate 
would still be true for a given to, we would nevertheless lose uniformity of the constants in the 
estimates with respect to to- We already stress at this stage that it is this uniformity that will 
be crucial in the proof of Proposition [ 



Proof of Lemma \3.6l It will be convenient to work with the function 



w := t x w 
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for some smooth periodic function A : U — > R, in order to optimize the positivity requirement on 
the matrix N at the end of this proof. Since w is a smooth solution of (I3.8[) . 

D 2 w + 2aDw + bw- t 2 k 2 d 2 x w = F[w] 

for t > to with coefficients a, 6, (or Ai and A2), and k, it follows by direct computation that 

D 2 w + 2aDw + bw- t 2 k 2 d 2 x w = t x f + F L [w}. 

Here, F^lw] is an expression linear in r/w, Dw and tkd x w with smooth coefficient functions, 
which arqj 0(< M ) at t = for some constant // > 0. Hence, the new source-term is again of the 
form (|3.9|) with the conditions (|3.10p . The coefficients of the principle part are given by 



De[w] = rj 2 wDw + DwD 2 w + (1 + Dk/k)(tkd x w) 2 + t 2 k 2 d x wd x Dw 
= rj 2 wDw + Dw{-2aDw ~bw + t 2 k 2 d 2 x w + t x f + F L [w}) 

+ (1 + Dk/k)(tkd x w) 2 + t 2 k 2 d x wd x Dw 
= -(b-r/ 2 )wDw - 2a(Dw) 2 + (1 + Dk / k){tkd x w) 2 + Dw{t x f + F L [w\) 

+ t 2 k 2 d 2 x wDw + t 2 k 2 d x wd x Dw. 

When De[w] is multiplied with t 2 ^ 2 ^) - "), the last two terms can be treated as follows 



The first term on the right vanishes after integration in space by virtue of periodicity on the 
domain U. Now we collect all terms of D E as follows 




DE[w] = - KTfb 1 E[w\ + e- Kt " f 2(3?A 2 - a)t 2 ^- a ^ <e[w]d: 




Now, 



t 2 ^ x ^- a+1 \k 2 d 2 x wDw + k 2 d x wd x Dw) 

= d x {t 2 ^ j2 ^- a+1 h 2 d x wDw) - 2t 2 ^ x ^ a \td x k){tkd x w)Dw 
-2t 2 ^ 2 ^- a \d x RX 2 ){tklnt)(tkd x w)Dw. 



DE[w) -. DEi [w] + DE 2 [w], 



where 



DE x [w\ : 




In the case that A is not a constant, this is strictly speaking only true if A is not too negative. 
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DE 2 [w] := ie" K * T 



and 

t 2{K\ 2 - a ) ^ _ K7t 7(^)2 _ Klt 1( D fif _ K-ftf (tkd x w) 2 

+ 2(t x f Q + F L [w])Di5jdx. 

Using the expressions of Ai and A2, we get 

DEi[w\ := J f 2 (^-A-a) ^ _ (a - 5RA 2 + a^wf - (RA X — A + a)(Dw) 2 

Dk 

- (A - 3?A 2 + a - 1 - — )(tkd x w) 2 

k 

- ((3Ai) 2 /t7 — 7/) (i]w)Dw 

- 2(ta E fc + 9 X (3?A 2 - A) (tfc In t))(tkd x w)Dw) dx. 

When we choose A = 5RAi, as we will do now, we can write DE\ as follows 

DE 1 [w\= / -(%■ N -w 1 )dx, 
Ju 

where N is the energy dissipation matrix in (|3.1ip and 

~ . = f SRA 2 (x)-a^~^ ) x^DQfa x ) t tk(t, X)d x w(t, x)). 

Note that w — w ■ T with w from (|3. 14|) and 

/l SRAi/t? (a a .«A 1 )(tfclnt)M 
T := 1 . (3.17) 

\o 1 y 

The matrix T is invertible and can hence be interpreted as the transformation matrix from the 
variable w to the variable w and vice versa. Thus, if TV is positive semidefmite at all (t, x), it 
follows that DE\\w\ < for all w, and hence for all w. According to the hypothesis of this 
lemma, this is the case, and we are left with 

DE[w](t) < DE 2 [w](t). 

Denote by (•, -)l 2 (u) the Euclidean L 2 -scalar product. Thanks to the properties of the expression 
Fl[w], we can choose n > (large enough) and 7 > (small enough) so that, uniformly for all 
to (provided 5 is sufficiently small), 

DE 2 [v^(t) <e- Kt " {(0,t B ^- a f o (t,-),0)M,-))L H u). 

In order to integrate the differential inequality for E[w] now, we use p. 151) and the Cauchy- 
Schwarz inequality to obtain 

D(^e-^\^(t,-)\\l 2(u ^ <e-^\\t^-y a (t r )\\ L2{u) \\^(t,-)\\ LHU) . 

This yields 

j t \\w(t, -)\\ LHU) < i/^^IRt, -)\\ LHU) + r 1 ||t RAa - '/o(* J OIUw 
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Then the Gronwall inequality implies 

for all to < t < 8. Thanks to the properties of the transformation matrix T defined in (|3.17D 
under our assumptions, one can check that there exists a constant C > (independent of to and 
t) so that 

\\G(t, OlU^a) < CeW-W (\\w(t , -)\\» (u) + J* 8- l \\a^'-" f (s, ^h^ds 

□ 

Proof of Proposition \3.4\ We start by assuming a smooth function /o on t > 0, and look for 
weak solutions w £ Xs. a ,i- The first step is to consider a monotonically decreasing sequence 
(Tn)neN C (0,5] converging to 0. We define a sequence (w n )n£H of functions on (0,8] x U as 
follows. For all n £ N, we set w n (t) — for all i G (0,r„]. On the interval [r n ,(5], we let u>„ be 
the unique solution of the RIVP for zero regular data at to = r n . The linearity of the equation 
and the conditions on the coefficients imply that w n is well-defined on the whole interval (0,(5], 
and that w n g C 1 ([0,<5] x U). It is easy to see that indeed, (w n ) C Xs, a .i for all a > 0. The 
motivation for choosing the sequence (w n ) is that the associated functions v n (defined according 
to (|3.6p ) can be hoped to behave more and more like a solution to the second-order Fuchsian 
equation obeying the leading-order behavior dictated by u, when n tends to infinity; recall that 
our assumption u = represents no loss of generality. Hence, the sequence (w n ) is expected to 
converge to a solution w of the SIVP. We prove that this is the case making use of the energy 
estimates for the RIVP derived in Lemma \WM 

Fix some arbitrary m, n £ N with m > n, thus < T m < r n < 8, and define £ := w m — w n . 
Thus, £ is identically zero on (0, r m ], it satisfies (|3.8[) with the given source-term for [r m , r„], and 
it satisfies (|3.8p with the given source-term but with vanishing inhomogeneity /o for [r n , 8] (due 
to the linearity of the equation). Then, Lemma 13.61 implies that 

f=0, t£(0,T m ], 
U(t,-)\\mu) I <Ce^-^)f T s-^\\s^-<*h(s,-)\\ L * {u) ds, te[r m ,r n ], (3.18) 

(<Ce^-^Q iS - 1 \\ s ^- a f ( s ,-)\\ LHu) d s , teM, 

where £ is the vector- valued function associated with £ in the same way as in (|3.14p . In particular, 
all constants here are independent of r m and r n . We get 

s- e \\s^- a f (s, -)\\ L 2 (U) < ||/o|ka +e ,0, (3.19) 

for all s £ (0, 8]. Here, e > is the constant given in the hypothesis of this proposition. In total, 
the map s M> s _e ||s SR:A2_a /o(s, ■)IIl 2 (c/) i s a bounded continuous function on (0,8] and hence is 
integrable. Hence, the function 

G(t) := [ s- l \\s^- a f ( S ,-)\\ LHu) ds 
Jo 

is well-defined and finite for all t. Indeed, it is continuous for t £ (0, 8] and lim f ^o G(t) — 0. This 
implies that G is uniformly continuous on (0,(5]. By taking the supremum in t on the interval 
(0,(5] of p,18[) and adapting the constant C if necessary, we show that 

\\wm - WnWlo.,1 < C\G(r n ) - G(r m )|. (3.20) 
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In particular, we point out that while C can depend on the choice of 5, it is bounded for small 5. 
Now, since G is uniformly continuous and (r n ) is a Cauchy sequence with limit zero, it follows 
that (G(r„)) is a Cauchy sequence with limit zero. Thus, p. 201) implies that the sequence (w n ) 
is a Cauchy sequence in (Xg <a> ±, |j • ||^ Q x ), and hence there exists a limit function w S Xs. a ,i- 

Now we have to check that w is a weak solution of the Fuchsian equation, i.e. (V [w] , <p) = 
for all test functions (f> according to Q3.7[) . Consider an arbitrary test- function <j) and pick up any 
n € N. The sequence element w n is constructed so that 

ds. 



{V[w n U)\< ^ (s^- a f (s,-),<f>(s,-)) L2iu) 



This estimate holds since, for all t > r n , the approximate solution u>„ satisfies the equation 
This yields 

-1 II „»A 2 -a 

'"|£»(17)' 



|(7>„U)|< sup W{t,-)\\ LHU) / S - 1 || S ^- a /o(s,-)|| i3fm rfs = CG(r n ) ) 
te(o,<5] Jo 



for some constant C. Hence, lim n _ ! . 00 (V[w n },(f>) = 0. Since (V[-],4>) is continuous on Xs, a ,i with 
respect to the norm || ■ \\£ al for any given test function (f> as noted earlier, it follows that 

(V[w],(p) =0. 

Hence w is a solution of the SIVP. 

Let us check that the limit w of the sequence (w n ) does not depend on the choice of sequence 
(t„), i.e. that our solution procedure yields a unique solution (which, however, is not guaranteed 
to be the only solution of the SIVP at this stage of the proof). Let (r„) be a monotonically 
decreasing sequence in (0, 8] with limit 0, (w n ) the corresponding sequence of approximate so- 
lutions and w its limit; the same for another monotonically decreasing sequence (r„) in (0,5] 
with limit 0, sequence of approximate solutions (w n ) and limit w. Now, we take the union of 
the two sequences (r„) and (r„) and sort the new sequence (t„), so that it becomes a monoton- 
ically decreasing Cauchy sequence with limit 0. Then (|3.20p shows that \\w n — w m \\s' a i — > for 
n, m — > oo. Hence 

||W - W\\s, a ,l < \\W - U; n ||^ a ,l + \\ w n - Wm\\s,a,l + \\ w ™ ~ ^m||^ a ,l ~> °» 

and so w — w. 

So far, we have shown that for all /o € Xs. a + e .o H C°° ((0, 5] x U) , there exists a weak solution 
u? € Xs. a ,i, and that u; is independent of the choice of sequence (r n ). Then the solution operator 

H : Xs, a+e ,o n C°°((0, 5] x (7) -»■ X s , a>1 , f ^ «;, 
is well-defined. It is clearly linear, and we derive the estimate (I3.12[) now. From (|3.20p . we get 

\H\7, a ,i<\M\7^i + CG(S). 

We can estimate || xui ||<s,q;,i as follows. Because Wi is a solution of the RIVP with zero regular 
data at to = t% , estimate (13.16[) yields 

\\Mt,-)\\mu) <Ce^-^ [ 8- 1 \\J a *- a Ma,-)\\ L , (V) d8 
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for all t G [ti,5]. For all t G (0,ri], we have ||u?i(t, -)lli/ 2 (t/) = 0- This shows that ||iui||,s,a,i < 
CG(S) with some adapted C, and hence, absorbing the factor 2 into the constant, we get 

IMIw < CG(S). (3.21) 

Now, the estimate f|3.19[) gives 

G(S) < ^1/o||£ a+ e,o- 

Using this together with (|3T2"T) yields (l3~T2")l on the subset Xs, a+e ,o n C°°((0, 5} x [/) of X S)Ot+et0 . 
Now we proceed in the same way as in the proof of Proposition 12.61 in order to extend this 
operator and the validity of (I3.12j) to the full space. For any given test function </>, the expression 
(V[w],(j)) is continuous with respect to w and fo (in their respective norms). This is sufficient to 
show that the extended operator H maps to a weak solution w. 

Finally, estimate p. 121) allows us to show uniqueness. Assume that we have two solutions w 
and w for the same source-term. Then w — w is a solution of the same equation with fo = 0. 
Hence, (|3.12p implies that 

and so uniqueness is established. □ 
3.4 Non-linear theory in the spaces Xg >a! 2 and X^ at00 

The general non-linear theory The well-posedness theory of the previous section, where 
we restrict to the space Xg l0t: i, has certain limitations. First, the statement that the solution 
of the Fuchsian equation w is an element of Xs^ a ,i yields particularly weak information about 
the behavior of the first spatial derivative at t — 0. It would be advantageous if we were able to 
prove the stronger statement w G Xs ya ,i, possibly under stronger assumptions. Second, it turns 
out that we need to require a Lipschitz property of the source-term for the general non-linear 
case which rules out natural non-linearities, for instance quadratic ones, if we only control the 
first derivatives of the solution. In the case of one spatial dimension, as we always assume in 
the whole paper, it is sufficient to increase regularity to the space Xg, a ,2- It is then clear how to 
proceed to Xs, a ,k with arbitrary k € N. Nevertheless, in some applications [7J, the space Xg, a ,k 
imposes too strong a restriction, due to the weak control of the highest spatial derivative. This 
problem can be avoided by formulating the theory in the space Xs, a ,<x- 

Lemma 3.7 (Existence of solutions of the linear singular initial value problem in Xg^p). Let 
us make the same assumptions as listed in the beginning of Section \3.°2\ with f a = for a = 1, 2, 3 
(for simplicity) . Then there exists a unique solution w G Xg. a ,2 of the singular initial value 
problem for given 5, a > provided: 

1. The energy dissipation matrix (|3.1ip is positive definite at each (t,x) G (0,<5) x U for a 
constant r\ > 0. 

2. The source-term function fo is in Xs^ a +e,i for some e > 0. 
Then, the solution operator 

H : X s , a +c,\ X s . a ,2: fo H' w, 
is continuous and there exists a finite constant C e > so that 

P[/o]||£«, 2 <m||/olk Q+£ ,i, 
for all fo. The constant C e is bounded for all small 5. 
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Analogous results hold for systems and for general linear source terms of the form (|3.9[) with 
non- vanishing functions /i, fi and fa obeying decay conditions analogous to (13.101) also for the 
first derivatives. Moreover, the result can be generalized to an arbitrary number k of derivatives, 
i.e. to solutions in the space Xs^ a ,k- For k > 3 in one spatial dimension, the Sobolev inequalities 
imply that the weak derivatives can be identified with classical derivatives. Hence the solution 
w G Xs. ay 3 of the weak form of the equation fl3.7[) is then a classical solution of (|3.1[) with v = w, 
i.e. u = 0. 

Proof. One sees immediately that (the generalization to systems of) Proposition 13.41 applies 
directly to the system of equations for the unknowns (it>o, u>i, W2) with wo := w, w\ — Dw and 
W2 := t ei d x w, where t\ > can be any sufficiently small constant. We cannot choose e\ = since 
this would lead to a source-term which is not consistent with the hypotheses of Proposition [33] 
in the following. For this system the energy dissipation matrix has the required properties, if it 
has the required properties for the original equation for w and if we assume the same constant 
a for the equations for wq, w\ and W2 • However, the energy dissipation matrix must positive 
definite instead of positive semidefinite due to the presence of the non- vanishing constant t\. 
One obtains existence and uniqueness in a space Xs, a ,i for vector-valued functions (u>o, u>i, u^)- 
The thus obtained space Xg, a ,i for vector-valued functions (uio,wi,W2) equals the space Xs,Si,2 
for the original scalar function w where a differs from a by the arbitrarily small constant e%. □ 

It is important to note that, when we repeat the proof for an arbitrary number of derivatives 
k, the quantity a can be chosen arbitrarily close to a irrespective of the choice of k. 

Proposition 3.8 (Existence of solutions of the non-linear singular initial value problem in 
Xs >a ,2)- Suppose that we can choose a > so that the energy dissipation matrix (|3.11[) is positive 
definite at each (t, x) G (0, S) x U for a constant n > 0. Suppose that u = and that the operator 
F has the following Lipschitz continuity property: For a constant e > and all sufficiently small 
5, the operator F maps X$, a .2 into Xg <a + e i and, moreover, for each r > there exists C > 
( independent of S) so that 

\\F[w] ~ F[w]\\ s , a+etl <C\\w- iB\\s, a ,2 (3.22) 

for all W, W € B r (0) C Xg^ a ^- Then, there exists a unique solution w € Xs, a ,2 °f the singular 
initial value problem. 

The generalization of this result to arbitrarily many derivatives is again straightforward. 

Proof of Provosition \3.8[ Similar to Section [21 we define the operator G := H o F, and argue in 
the same way as in Proposition 12.71 and Theorem 12.81 that under the hypothesis, this operator 
is a contraction on closed and bounded subsets of Xs, a ,2 if 5 is a sufficiently small. Hence the 
iteration sequence defined by lOj+i = for j >1 and, say, w\ = converges to a fixed point 

w G Xs j0 i,2 with respect to the norm || • ||^ Q 2 . Because of the properties of H, a fixed point of G 
is a solution of the SIVP. Hence, we have shown existence of solutions. Uniqueness can be shown 
as follows. Given any other solution w in Xs^ a .2, it is a fixed point of the iteration Wj+i — G[wj]. 
Because G is a contraction, there, however, only exists one fixed point, and hence w = w. □ 

Proposition 3.9 (Existence of solutions of the non-linear singular initial value problem in 
Xs,a,oo)- Suppose that we can choose a > so that the energy dissipation matrix p. lip is 
positive definite at each (t,x) G (0,(5) x U for a constant n > 0. Suppose that u = and that 
the operator F has the following Lipschitz continuity property: For a constant e > 0, every 
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sufficiently small 8 > and every non-negative integer k, the operator F maps Xs t0!t k+i into 
Xs, a +e,k and, moreover, for each r > 0, there exists C > (independent of 8) so that 

\\F[w] - F[w]\\ SM <C\\w- w\\Z Q:k+1 (3.23) 

for all w,w G B r (0) n Xs, a ,k+i C Xs, a ,k+i- Then, there exists a unique solution w G Xs, a ,oo of 
the singular initial value problem. 

Here, B r (0) is defined with respect to the norm || • || £ a k+1 . We note that the constant C is 
allowed to depend on k. Note that the Lipschitz estimate involves the norm || • k+1 , while the 
elements for which this estimates needs to be satisfied are required to be only in the subspace 
Xs. a ,k+i of Xg, a .k+i- The main advantage of this result over the finite differentiability case is 
that we only need to check that F maps Xs^ a ,k+i into Xs^ a +e,k for all fc, instead of the stronger 
statement that F maps Xs, a ,k+i into Xs_ a +e.k (which would of course, however, only need to 
hold for finitely many k). 

Proof. We first generalize Lemma [3.71 to the case fo G Xg M:00 . Then the solution of the linear 
equation H[/o] G Xs, a .oo- Since F maps Xs^ a ,oo to itself, the operator G := HoF maps Xs, a ,oo to 
itself. Hence, the same iteration as in the proof of Proposition 13.81 leads to an iteration sequence 
(vjj) C Xs_ a ,oo- Let k > 1 be arbitrary. According to the hypothesis, F maps Xs ia ,k into 
Xs, a +e,k-i- Lemma 13.71 implies that H maps Xs^ a +e,k-i to Xs, a ,k' Hence G can be consider as 
a map Xs, a ,k — > Xg t a,k- We can consider X$ l0lt k to be a subset of Xg. a ,k- Then, if we choose 
8 small enough, the restriction of G to Xs, a ,k H B r (0) C Xs,a.,k is a contraction with respect 
to the norm || • ||j Q k due to (|3.23p . This implies that for an appropriate choice of r, we have 
that (uij) is a Cauchy sequence in Xs t a,oo H B r (0) C Xs, a ,k with respect to the norm || ■ W^ 1 k . 
Hence, it converges to a limit wr^ G Xg( k ^ a k . We have written 8(k) now instead of 8 in order 
to stress that 5 does depend on k. We get such a limit function w^) G X$/f.),a,k f° r a U k > 1. 
It is straightforward to check that wr^ty) — w (k 2 )(t) for all two integers k\,k2 > 1 for all 
t G (0, min{<5(fci), (5(^2)}]. Hence, it follows that wqa G Xstk+i),a,k (without tilde!) for all k, 
and this means that wr^) € Xs(k), a ,k after possibly decreasing 8(k) > sufficiently for all k. So 
for any given k, the limit W(h\ is in the range of the operator G, and thus u)\k) 1S ^ ne unique 
fixed point of G and so the unique solution of the equation in Xg( k j a k . However, it is not 
obvious at this point whether we are forced to choose 8(k) — > for k — > oo, and we are left 
with demonstrating that this is not the case. As soon as we have this, we have constructed 
the solution w G Xg^ attyo for some 8 > 0. This, however, requires only standard arguments for 
symmetric hyperbolic equations. Hence, we find that we can choose 8 — 8(3) (for one spatial 
dimension) . □ 



The (standard) singular initial value problem The following discussion is devoted to 
particular choices of the function u motivated by the heuristics introduced in Section [2] Consider 
the case that u is given by 

fu*(x)t- a Wlnt + u**(x)t- a W a 2 = b, 
u(t,x) —un(t,x) :— < . , , . , , „ (3.24 

V ' K ' \u*(x)t- Xl(x) +u.*{x)t- x »( x \ a 2 ^b, 1 ' 

with u Sf ,u*i t given by (|2.7p and with asymptotic data G H 3 (U). In this case, we will 

speak of the standard singular initial value problerrQ Note that this means that for all 
t > 0, the map u(t, ■) and all its time derivatives are in H 3 (U). 



In order to simplify the language, we often speak of the singular initial value problem if there is no risk of 
confusion. 
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Theorem 3.10 (Well-posedness of the standard singular initial value problem in Xs,a.2)- Given 
arbitrary asymptotic data u*,u** G H 3 (U), the standard singular initial value problem admits 
a unique solution w £ Xs, a .2 for a, 5 > 0, provided S is sufficiently small and the following 
conditions hold: 

1. Positivity condition. Suppose that we can choose a > so that the energy dissipation 
matrix (|3.11l) is positive definite at each (t,x) € (0,(5) x U for a constant r\ > 0. 

2. Lipschitz continuity property. For the given a > 0, the operator F satisfies the Lipschitz 
continuity property stated in Provosition \3.8\ for all asymptotic data w*,u*» G H 3 (U) for 
some e > 0. 

3. Integrability condition. The constants a and e satisfy 

a + e<2(^(x) + l)-SR(Ai(x)-A 2 (x)), x € U. (3.25) 

An analogous theorem can be formulated for the C°°-case based on Proposition 13.91 In 
this case, the asymptotic data it*, u** must be in C°°(U) and the Lipschitz condition must be 
substituted by the condition of Proposition 13.91 The unique solution w of the singular initial 
value problem is then an element of Xs t a,oo- 

We note that there might be room for improvements in the finite differentiability case k = 2, 
since three derivatives of the asymptotic data yield control of only two derivatives of the solution. 

Proof. We can apply Proposition 13.81 if we are able to control the additional contribution of the 
term L[u] which has to be considered as part of the source-term. It has no contribution to the 
Lipschitz estimate (|3.22[) , but we have to guarantee that under these hypotheses, L[u] G Xs, a +e,i 
for the given constant e. This is indeed the case if (|3.25p holds. □ 

Example 3.11. Consider the second-order hyperbolic Fuchsian equation 

D 2 v - XDv - t 2 d 2 x v = 0, 
with a constant A. This is the Euler-Poisson-Darboux equation. In the standard notation it is 

d 2 t v-d 2 x v = - t {\-l)d t v. 

Note that \ = 1 is the standard wave equation, and in this case, the standard singular initial 
value problem reduces to the standard Cauchy problem. 

1. Case A > 0. With our notation, we have \\ = 0, A 2 = — A, /3 = 0, v = 1 and f = 0. 
The positivity condition of the energy dissipation matrix p. lip is satisfied precisely for 
a > 1 — A and all sufficiently small rj > 0. The integrability condition (13.251) is satisfied 
precisely for A < 2 — a. Hence, our previous proposition implies that the singular initial 
value problem is well-posed, provided 

< A < 2. (3.26) 

Namely, in this case there exists a solution w in Xs, a ,2 for some a > for arbitrary 
asymptotic data in H 3 (U). 
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2. Case A < 0. With our notation, we have Ai = |A|, A2 = 0, j3 = 0, v = 1 and f = 0. 
The positivity condition of the energy dissipation matrix (|3.11[) is satisfied precisely for 
a > 1 — |A| and sufficiently small 7/ > 0. T7te integrability condition (I3.25[) is satisfied 
precisely for |A| < 2 — a. Hence, our previous proposition implies that the singular initial 
value problem is well-posed, provided 

-2 < A < 0. 

Namely, in this case there exists a solution w in Xg, a ,2 for some a > for arbitrary 
asymptotic data in H 3 (U). 

Now, it turns out that general smooth solutions to the Euler-Poisson-Darboux equation can be 
expressed explicitly by a Fourier ansatz in x and by Bessel functions in t. It is then easy to 
check that (|3.26[) (and similarly for A < 0) is sharp: While for < A < 2, all solutions of the 
equation behave consistently with the two-term expansion at t = 0, this is not the case for A > 2 
for general asymptotic data. Hence the singular initial value problem is not well-posed for A > 2. 
This is completely consistent with our heuristic discussion in Section \2.3\ Namely, if X = 2, the 
assumption that the source-term t 2 d 2 v is negligible a£ i = fails since it is of the same order in t 
att~0 as the second leading-order term. However, we can see in the proof of Theorem V3.HA that 
in the special case it* = (and arbitrary u**), the integrability condition (|3.25p can be relaxed. 
For this special choice of data, solutions to the singular initial value problem exist even for A > 2. 

Singular singular initial value problems with asymptotic solutions of order j One 

of the main aims of this paper is to study the well-posedness of the standard singular initial 
value problem just discussed. In this sense, we can be satisfied with Theorem 13.101 However, 
it turns out that, often in applications, the three conditions in this theorem cannot be satisfied 
simultaneously. While it is often possible to find constants a and e in accordance with the second 
and third condition, it can turn out that the corresponding choice of a is too small to make the 
energy dissipation matrix positive definite. The following trick can sometimes solve this problem. 

For the following discussion, we need to bring together results from Sections [5] and [3J and we 
are forced to distinguish between operators now which for the sake of simplicity have carried the 
same name so far. Consider some asymptotic data it* and it** and define the function ito as in 
p. 241) . We write / for the source-term in (|2.ip and continue to write / for the source-term in 
(|3.1j) . Accordingly, we write F[w] := f[uo + w] and F[w] := f[uo + w], so that e.g. 

F[w] = F[w] + t 2 k 2 d 2 x {u + w). 

In the same way, we distinguish between operators L and L. Now we make the same assumptions 
on F as in Proposition 12 . 71 in Section l2~5l If the asymptotic data is in H mi (U) for some positive 
integer mi, then the function Uj, referred to as uij in Proposition 12 .71 with iti — is well-defined 
in Xs,a,i,mx-2{j-X) f° r some a > and all j with mi/2 + 1 > j > 1, provided F maps Xg^a,i,m 
to Xs,a+e,l-i,m-i f° r all m < mi, for some integer I > 1 and for some e > 0. 
Now let us choose the leading-order function it as 

u(t, x) = Uo(t, x) + Uj(t, x), (3.27) 

for j in the range given above. We refer to the singular initial value problem based on this 
choice of leading-order term as singular initial value problem with asymptotic solutions 
of order j. For j = 1, it reduces to the standard singular initial value problem; hence we will 
focus on the case j > 2 in the following. Note that, if w is a solution of the singular initial value 
problem of order j, it is also a solution of the standard singular initial value problem. However, 
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if there is only one solution w of the singular initial value problem with asymptotic solutions of 
order j for given asymptotic data, it does not mean that w is the only solution of the standard 
initial value problem for the same asymptotic data. 

It can be seen easily that the remainder w of a solution of the singular initial value problem 
with asymptotic solutions of order j satisfies the equation 

L[w] = Fj[w] := F[uj + w] - F[uj-x] + t 2 k 2 dl{uj - 

for j > 2. Now thanks to Theorem 12. 101 we have 

U j ~ U j-1 ^ ^S,a+{j-2)Ke ,l, mi -2{j-l)j 

for all k < 1. Hence it is reasonable to restrict to remainders w £ -^5,a+(j-2)«eo,;,mi— 2(j-i)— l 
in the following. One finds easily that this means that Fj[w] £ Xsa+(j-i)Ke ,l-i,mi-2(j-i)-2 if 
2(j3(x) + 1) > ne for all x £ U by using similar arguments as in Theorem 12.101 This gives us 
hope that we can apply Proposition 13.81 with 

a := a + (J - 2)«e. (3.28) 

The effect of the ansatz (|3.27[) is a value of a which increases a by (j — 2)kc. Namely if mi is 
sufficiently large, we can choose j large enough so that the energy dissipation matrix, evaluated 
for a, can become positive definite. The main prize that we pay with this approach is that the 
asymptotic data must be sufficiently regular and that we must live with a loss of regularity which 
is stronger the larger j is. 

For the statement of the following theorem, we need the following notation. For all w £ Xs, aj k 
(or w £ Xg.a.k respectively), we introduce the functions Es, a ,k[w] ■ (0,(5] — > K (or Es ia .k[ w ] '■ 
(0,(5] — >• R respectively) which arc defined in the same way as the respective norms, but the 
supremum in t has not been evaluated yet. In particular, this means that Eg iaj k[w] (or E$ tCXj k[w\) 
is a bounded continuous function on (0,(5]. 

Theorem 3.12 (Well-posedness of the singular initial value problem with asymptotic solutions 
of higher-order in Xs ia ,2)- Given any integer j > 2 and any asymptotic data u*,?/*, £ H mi (U) 
with mi = 2j + 1, there exists a unique solution w £ Xs, a .2 of the singular initial value problem 
with asymptotic solutions of order j for some a > 0, provided 

1. F maps Xs,u,mx into Xs,a+e,mx-i f or a H asymptotic data u*, u** £ H mi (U) for some e > 
and a given by (|3. 281) for an arbitrary k < 1 . 

2. The characteristic speed satisfies 

2(j3(x) + 1) > ne for all x£~U 
for the same constant k chosen earlier. 

3. F satisfies the following Lipschitz condition: for each r > there exists a constant C > 
( independent of 5) so that 

E s ,a+e,i[F[w} - F[w]](t) < CE s ^ 2 [w - w](t) 

for all t £ (0, 6] and for all w,w £ B r (0) C Xs,a. 2 • 

4- The energy dissipation matrix (13.111) (evaluated with a) is positive definite at each (t,x) £ 
(0, S) x U for a constant n > 0. 
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The third condition above is meaningful since both sides of the inequality are continuous and 
bounded functions on (0,6]. Note that this theorem can be formulated without difficulty for the 
C°°-case and leads to a simpler statement. 

Proof. Under the first hypothesis, it can be shown similar to Proposition ^ . 7l that the j-th element 
of the iteration sequence Uj is in Xg & mi _ 2 (j-i)- As in Theorem 12.101 it can be demonstrated 
that u - Uj-i is in X s ,a+(j-2) Ki , mi -2(j-i) for all j > 2. Now let w e X g ,a+(j-2) K e,mi-2(j-i)-l- 
Under the second hypothesis, it follows that -Fj[W] G A ( 5 l 5+j Ke .mi-2(j-i)-2- Now, in order to 
apply Proposition 13.81 with F substituted by Fj, it is necessary to choose mi = 2j + 1. The 
operator Fj satisfies the Lipschitz condition of Proposition [3TH] if the third hypothesis is satisfied. 
Thanks to the fourth assumption, we can now apply Proposition 13.81 □ 
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